Two distinct Mott-Insulator to Bose-glass transitions and breakdown of self averaging 

in the disordered Bose-Hubbard model 
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We investigate tlie instabilities of the Mott-insulating phase of the weakly disordered Bose- 
Hubbard model within a renormalization group analysis of the replica field theory obtained by 
a strong-coupling expansion around the atomic limit. We identify a new order parameter and as- 
sociated correlation length scale that is capable of capturing the transition from a state with zero 
compressibility, the Mott insulator, to one in which the compressibility is finite, the Bose glass. 
The order parameter is the relative variance of the disorder-induced mass distribution. In the Mott 
insulator, the relative variance renormalizes to zero, whereas it diverges in the Bose glass. The di- 
vergence of the relative variance signals the breakdown of self-averaging. The length scale governing 
the breakdown of self-averaging is the distance between rare regions. This length scale is finite in 
the Bose glass but diverges at the transition to the Mott insulator with an exponent of = 1/D for 
incommensurate fillings. Likewise, the compressibility vanishes with an exponent of7 = 4/_D — 1 
at the transition. At commensurate fillings, the transition is controlled by a different fixed point at 
which both the disorder and interaction vertices are relevant. 

PACS numbers: 05.30.Jp, 72.15.Rn, 64.70. Tg, 64.60. ae 



I. INTRODUCTION 

Two distinct mechanisms localize interacting bosons 
moving in a random environment. Either repulsions from 
local on-site interactions on a D-dimensional lattice lock 
the bosons in place or disorder inhibits tunneling from 
site to site. The former insulating state is an incompress- 
ible Mott insulator (MI) , while the second is a compress- 
ible Bose glass (BG). In this paper, we identify an order 
parameter and an associated length scale that can de- 
scribe the vanishing of the compressibility at the MI/BG 
transition. The order parameter is the relative variance 
of the disorder-induced mass distribution. The key in- 
sight leading to this observation is that despite the rele- 
vance of the disorder vertex in D < 4 spatial dimensions, 
it is possible to identify a Mott-insulating fixed point 
where the mass, corresponding to the Mott gap, diverges 
fast enough such that the relative variance of the mass 
distribution is renormalized to zero. Sufficiently close to 
the MI, both the mean and the relative variance diverge, 
identifying the adjacent phase as an insulating BG, irre- 
spective of the boson filling. We will show that the length 
scale underlying this order parameter is the effective dis- 
tance between rare regions. 

In their initial treatment of this problem, known as 
the disordered Bose-Hubbard (BH) model. Fisher and 
colleagues^ situated the BG between the MI and super- 
ffuid (SF), thereby preventing a direct SF/MI phase tran- 
sition. They pointed out however, that it is in principle 
possible but extremely unlikely that for sufffciently weak 
disorder the BG phase is completely suppressed at com- 
mensurate boson fillings. For more than twenty years, 
this question remained controversial as simulations and 
analytical arguments both support2r2i and negatei^"— 
the possibility of a direct MI/SF transition in the pres- 



ence of weak disorder. Only recently, strong arguments^^ 
based on a mathematical theorem have been presented 
for the existence of the BG upon the destruction of 
the MI, thereby precluding a direct transition to the 
SF. Given this result, it is important to establish pre- 
cisely how the MI to BG transition obtains. In this pa- 
per, we are concerned with the critical theory for the 
MI/BG transition for strongly interacting bosons on a 
D-dimensional lattice subject to weak disorder. Of par- 
ticular interest is the universality of the transition, e.g. 
the identification of a diverging length scale and the de- 
termination of the critical exponent describing the diver- 
gence of the inverse compressibility. Despite extensive 
numerical2i^i^i^i^i^ii^ and analytical^ii2,-i£ii^-^ studies 
of the disordered BH model, the nature of the MI/BG 
transition remains elusive. 

To be more specific, in the following we consider the 
disordered BH model in its simplest form. 



H = -i^(6|S, -|-h.c.)-f ^(e, 



(1) 



describing bosons with creation and annihilation opera- 
tors bl, bi (hi = blbi) which are hopping with amplitude t 
between neighboring sites i,j on a 15-dimensional hyper- 
cubic lattice and which are subject to an on-site repulsion 
U. The chemical potential /i controls the boson filling 
and disorder enters via random uncorrelated on-site po- 
tentials ei which must be appropriately bounded for the 
model to exhibit stable MI phases>i 

Since the formation of the MI is a consequence of 
strong repulsive interactions, the MI/SF transition in the 
clean system can only be understood within an effective 
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long-wavelength theory which is dual to the original BH 
model (dl and derived by a strong coupling expansion 
around the atomic limit.— >2ii^ The mass in the theory 
is given by the Mott gap which vanishes at the transi- 
tion and turns negative in the SF. In the latter broken- 
symmetry state the finite bosonic order parameter is pro- 
portional to the SF density. Although the underlying 
strong coupling expansion breaks down in the SF where 
the renormalization group (RG) flow is towards weak 
coupling,™ the effective long-wavelength theory can be 
used to analyze the instability of the MI. Since the theory 
is controlled in the MI state it should serve as a starting 
point to also analyze the instability towards the forma- 
tion of the BG in the presence of weak disorder. 

At the outset, it would seem that an effective theory 
in terms of the bosonic order parameter is completely 
inapplicable to the MI/BG transition since in both in- 
sulating phases the mass gap is finite and the SF corre- 
lations are short ranged. Certainly, the relevant length 
scale at the MI/BG transition has nothing to do with the 
SF correlation length but instead should be intrinsically 
related to the disorder. Such a divergent length scale 
should exist since the transition is expected to be con- 
tinuous with a divergent inverse compressibility. In this 
work, we demonstrate that the order parameter for the 
MI/BG transition can be extracted from the disorder- 
averaged replica version of the bosonic order parameter 
theory. The relevant quantity is the relative variance of 
the disorder-induced mass distribution which renormal- 
izes to zero in the MI and diverges in the BG and hence 
serves as the order parameter for the MI/BG transition. 
The corresponding length scale is the distance between 
rare regions in the BG which causes the breakdown of 
self-averaging in the system as indicated by a diverging 
relative variance. This correlation length is finite in the 
BG and diverges on approaching the MI. 

From the analysis of the RG fiow of the mean value and 
relative variance of the mass distribution we obtain the 
key results that (i) the instability of the MI is always to- 
wards a BG, in agreement with recent arguments,— and 
that (ii) the MI/BG transitions at incommensurate and 
commensurate boson fillings are not in the same univer- 
sality class. In the former case of incommensurate fill- 
ings, we calculate the critical exponents at 1-loop order: 
V = 1/D for the correlation length and 7 = 4/Z? — 1 for 
the inverse compressibility. 

The paper is organized as follows. In Sec. [TTl we argue 
that the relative variance of the disorder-induced mass 
distribution serves as the order parameter for the MI/BG 
transition which vanishes in the MI and acquires a finite 
value in the BG. We further identify the corresponding 
correlation length as the distance between rare regions in 
the BG which cause a breakdown of self-averaging in the 
system. The disorder-averaged replica theory describing 
the long- wavelength physics at strong coupling is derived 
in Sec, mil In Sec. lIVi the RG equations are reviewed and 
reformulated in a new set of variables which permits us 
to extract the scale dependence of the relative variance 



of the induced random-mass distribution. The results 
are presented in Sec. |Vl The phase boundary between 
the MI and the BG is determined by a numerical inte- 
gration of the RG equations. For incommensurate boson 
fillings, we obtain analytical results for the correlation 
length and compressibility exponents. Finally, we show 
that the MI/BG at commensurate boson fillings is in a 
different universality class. In Sec. I VI I our main results 
are discussed. 



II. BREAKDOWN OF SELF-AVERAGING 

The MI in the disordered BH model exists in the strong 
coupling limit. While the disorder must be appropriately 
bounded for the MI to persist,- the key limit that defines 
strong coupling is the ratio U/t ^ 1, where U is the on- 
site energy and t the hopping matrix element. As with 
all strong coupling problems, the natural variables that 
uncloak the physics are not related straightforwardly to 
those in the UV-complete Hamiltonian. Further, if a crit- 
ical theory is to correctly describe the destruction of the 
MI, it should contain the seeds of the BG phase, which 
Fisher and co-workers^ argued are due to the physics of 
rare regions. Consequently, the underlying critical theory 
might have nothing to do with the bare bosonic propaga- 
tor and the associated superfluid correlation length but 
rather a new length scale that is intrinsically related to 
the disorder. 

Consider the atomic limit of this problem. A scal- 
ing analysis around this regime^'^ demonstrates that the 
correlation length, ^, defined as the length scale beyond 
which the system encounters rare-region physics is finite 
in the BG and diverges at the transition to the MI with 
an exponent of = 1/Z? for generic bounded disorder 
distributions such as the box distribution P{e) — 1/(2A) 
for |e| < A and P{e) = for |e| > A. The violation 
of the lower bound of > 2/D known as the quantum 
Harris criterium?^ was attributed to a breakdown of self- 
aver aging2^ in the BG phase. That is, the fluctuations 
are not governed by the central-limit theorem. Several 
years ago, Aharony and Harris^^ showed that a break- 
down of self-averaging implies that the relative variance 
of any thermodynamic quantity must be finite. Taken in 
tandem, these results imply that characterizing a transi- 
tion governed by the physics of rare regions requires two 
quantities to be finite: 1) a finite length scale, ^, over 
which a rare region is encountered and 2) a finite value 
of the relative variance of any thermodynamic quantity. 

The central idea we advance here is that the nature 
of the MI/BG transition can be understood using per- 
turbative RG techniques but with a new set of variables 
and length scales that contain the information about the 
rare regions. The mass r in the effective long wavelength 
theory in the clean system corresponds to the Mott gap 
which renders the correlations of the bosonic order pa- 
rameter short ranged. Any form of disorder in the micro- 
scopic BH model, e.g. the potential disorder following 
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a distribution P{e), will induce an effective disorder dis- 
tribution P{r) of the mass coefficient r. Since in any in- 
sulating phase the superfluid order parameter must van- 
ish on large scales, in both phases, the MI and the BG, 
the mean r of the distribution must diverge under the 
RG, f{£) ^> oo. Hence, it cannot be used to distinguish 
between the MI and BG phases. Given the general con- 
siderations above, the natural quantity to distinguish the 
BG from the MI is the relative variance Rr of the mass 
distribution. 



Rr = (2) 

In the MI phase, the system is self averaging and there- 
fore Rr (i) since relative variances of extensive quan- 
tities have to vanish in the thermodynamic limit as a 
consequence of the central-limit theorem. Since the BG 
is characterized by a breakdown of self averaging, Rr is 
non-zero in the thermodynamic limit and therefore serves 
as the order parameter for the MI/BG transition. In this 
paper, break down of self-averaging refers entirely to the 
fact that Rr ^ 0. Note, Rr{l) ^ in the MI implies 
that ^ diverges at the transition; that is, there is no fi- 
nite length scale over which rare regions exist. We point 
out that a break-down of self averaging, as defined here, 
has nothing to do with a finite-size effect. Rather, it 
is a function entirely of whether or not in the thermo- 
dynamic limits ^, the distance over which rare regions 
exist, is finite. If it is, Rr is necessarily non-zero as we 
will demonstrate by an explicit calculation. 

To work out the RG flow of Rr-, we will use that the 
disorder vertex in the disorder-averaged replica theory 
is proportional to the variance of the mass distribution. 
We will proceed by recasting the RG equations in terms 
of Rr and show that these equations admit a stable MI 
fixed point at f = oo and i?^ = as well as a critical 
MI/BG fixed point at f = cxd and finite Rr. We point out 
that such an identification of fixed points is only possi- 
ble because of the identification of the relevant variables. 
Whereas the standard equations that are used to study 
this problem only capture a runaway flow away from the 
MI/SF fixed point of the clean system,A~ in the new 
variables we find that the unstable MI/SF fixed point is 
connected with the new critical MI/BG fixed point by 
a separatrix. Our analysis will further demonstrate that 
the MI/BG transition at commensurate fillings belongs 
to a different universality class. 

III. EFFECTIVE REPLICA THEORY 

We analyze the instability of the MI state towards the 
formation of the BG using the disorder-averaged replica 
version of the effective long-wavelength theory which is 
derived by an expansion around the atomic limit of the 
stable MI state. In this theory, the replica mixing dis- 
order vertex is proportional to the variance of the mass 
distribution which allows us to extract the RG fiow of 



the relative variance Rr ^ which we have identified as 
the order parameter for the MI/BG transition. In the 
absence of disorder, the effective action capturing the 
long- wavelength physics at strong coupling, tjl] ^ 1, 

\ (7irao0 + 72|ao<^|' + |V0|2-f r|0|2 + ft|0|4), 

(3) 

can be constructed using a coherent state path-integral 
representation of the BH model ([Ij and decoupling the 
off-diagonal hopping terms with a Hubbard-Stratonovich 
(HS) transformation. In this action, we have rescaled 
the length and imaginary time r to dimensionless units, 
X — Ar with A the momentum cutoff and xo = Vt. The 
precise dependence of the coupling constants on the mi- 
croscopic parameters is known analytically j22 For exam- 
ple, the mass coefficient which corresponds to the Mott 
gap, is given by 

/ m + 1 TO \ , , , 

r = 1 - y + — 4 

\TO — a; 1 — TO + xy 

with a; = /i/C/, 2/ = 2Dt/U ^ and ni the number of bosons 
per site. The interaction vertex h is proportional to re- 
fiecting the underlying strong-coupling expansion around 
the atomic limit (y = 0). An explicit expression for h 
can be found elsewhere. At the mean-field level, the 
MI/SF phase boundary is determined by r = 0. The MI 
states for different fillings m are characterized by a finite 
Mott gap, r > 0, and (0) = 0. The SF obtains for r < 
where the SF density is proportional to (</>) 7^ 0. 

As a consequence of the C/(l) symmetry of the fields?^ 
the temporal gradient terms are given by 71 = —dr/dx 
and 72 — —^d^r/{dx)^ and are therefore related to the 
slope and the curvature of the mean-field phase bound- 
ary. At the tips of the MI lobes at values Xm of the 
chemical potential where the boson filling is commensu- 
rate with the lattice and the MI states are most stable, 
7i = 0, and the critical theory is characterized by a dy- 
namical exponent z = 1. At incommensurate fillings, 
71 7^ 0, corresponding to z = 2 dynamics. 

Consider now the disordered case. The site energies 
will be uncorrelated and chosen from a bounded distribu- 
tion, £i g [—A, A]. The identical HS transformation can 
be performed^S yielding a dual theory with coefficients 
given by the expressions for the clean system but with 
the chemical potential x — fi/U shifted by the random 
potential ei/U on each site, e.g. = r{x — ei/U) for 
the mass term. This expression is well defined on every 
site if the disorder distribution P{ei) is bounded, and the 
chemical potential x lies in the intervals [m—1 + 6, m — S\ 
and 5 = A/U < 1/2. Note that this stability condition 
ensures that the system exhibits stable MI phases in the 
presence of disorder. 

The procedure for constructing the effective action is 
now standard: restore translational symmetry by per- 
forming the disorder average of the free energy using the 
replica trick and and take the continuum limit. In the 
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resultant expression, 



E 

a0 



(5) 

the first term corresponds to n identical copies of the 
disorder- averaged action with a the replica indices. 
Note that also for the disorder-averaged coefficients 71 = 
—dr/dx and 72 — ~^d'^f/{dx)'^. The replica-mixing dis- 
order vertex arises from the quadratic order of the cumu- 
lant expansion and is given by the variance of the mass 
distribution, 



g = {ArY = in - ff 



(6) 



As is evident, it is non-diagonal in imaginary time as 
a consequence of the perfect correlation of the disorder 
along the imaginary time direction. Since the disorder 
also couples to the temporal gradient terms, the replica 
theory contains disorder vertices with additional time 
derivatives.— These terms turn out to be irrelevant 
and do not change the universality of the transitions be- 
tween the MI and the BG and hence are not considered 
any further. 

As a consequence of the underlying strong-coupling ex- 
pansion, the variance (Ar)^ of the mass distribution is 
not simply given by the variance of the random-site dis- 
order in the BH model. This is easily seen by expand- 
ing Ti — r{x — €i/U) for small ei/U yielding (Ar)^ — 

-il'^/U^ + 7|(? - ^Vt^" + ■ • ■ with the n-th mo- 
ment of P{e) = P{—e). Hence, the effective disorder g 
is a continuous function of the chemical potential enter- 
ing via the coefhcients 71.2 and is strongly suppressed 
at the tips of the MI lobes where 71 = 0. Therefore, 
the inclusion of higher moments of the distribution P{t) 
is essential not only to account for the boundedness of 
the distribution but also to capture the coupling to dis- 
order at commensurate fillings. In order to retain the 
moments to infinite order, without loss of generality, we 
use the discrete distribution 

F(e) = (1 - 2p)5{€) + p5{e - A) + p5{€ + A), (7) 

corresponding to on-site energies which occur with the 
same probability p < 1/2 increased or decreased by A < 
U /2. The resulting mean r and variance g of the induced 
mass distribution entering the replica theory ([5]) are given 
by 



r{x-e,/U) 

(1 - 2p)r{x) +p[r{x -5)+ r{x - 5% 



(8) 



g = [r(a; -r)]2 

= (1 - 2py(x) +p[r'^{x - 5) + r2(a; + 5)] - r2,(9) 

and contain the moments of P(e) to infinite order. The 
MF phase boundary between the m = 1 MI and the SF 
obtained by the condition r = is shown in Fig. 1(a). 
Whereas at incommensurate fillings the interaction ver- 
tex h is irrelevant for _D > 2, the disorder vertex is always 



relevant in D < 4 rendering the MF theory meaningless. 
In order to determine the stability region of the MI and 
subsequent phase transitions, we employ the RG. 



IV. RENORMALIZATION GROUP ANALYSIS 

We proceed with the RG analysis of the effective 
replica field theory ([5]) . From the RG flow of the mean f 
and the variance g of the disorder-induced mass distribu- 
tion we can subsequently determine the scale dependence 
of the relative variance Rr = g /r^ to analyze the insta- 
bility of the MI towards the formation of a BG. 



A. Renormalization Group Equations 

After successively eliminating modes of highest energy 
corresponding to momenta from the infinitesimal shell 



-M 



< |fe| < 1 and rescaling of momenta (fc — > fee ), 



frequencies (w toe ), and fields 
obtain to one-loop order- 



-Ad£ 
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we 



dr 
d£ 
d7i 
d£ 
d72 
d£ 
dh 
d£ 
dg 
d£ 



2f+ 2/1/1 -Jog, 
(2 - z)7i + I^^ig, 
(2-2z)72 + /2(/o7?+72)5, 



(10a) 
(10b) 
(10c) 



(4 - D - z)h - (Jo + 4:Ji)h^ + QI^gh,{lOd) 
{A- D)g + Allf -AJifh, (lOe) 



where the scaling dimension A of the fields has been 
determined such that the coefficient of the spatial gradi- 
ent term remains constant. Further, we have absorbed 
a numerical factor arising from the shell integration by 
redefining h = 25i3/(27r)^7i and g = SD/{2n)^g with 
Sd the surface of the Z?-dimensional unit sphere. The 
frequency integrals are defined as Iq ~ C(0) = 1/(1 + r), 
h = LC(c.) = [7? + 472(1 + r-)]-V2, J, ^ lC^(u;) ^ 
/0/1/2, and Ji = I |C(^)|2 = 272/f with I = ^ 
and C{uj) = 1/(1 + r — i'jiuj + 72^^^) the on-shcU propa- 
gator. 



B. New Variables 

From the RG equations, it is evident that for D < A, 
the variance g increases with the scale which suggests 
that the RG fiow is always towards strong disorder. 
This interpretation is inconsistent with perturbative ar- 
guments which show that the disordered BH model ex- 
hibits stable MI phases in the regime of weak disorder,— 
and does not take into account that in the insulating 
phases the mean f also diverges under the RG. There- 
fore the question of whether disorder is strong or not 
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can only be answered by analyzing the relative variance 
Rr ^ g/^^ ■ In the stable MI phase, disorder is irrelevant 
and we expect that — 0. 

Realizing that f diverges for both insulating phases, we 
implement the approach of Aharony and Harris^- by fo- 
cusing on the relative variance i?,. of the disorder-induced 
mass distribution. As we argued in Sec. |TT1 serves as 
the order parameter for the MI/BG transition character- 
ized by a breakdown of self-averaging. 

Putting this into practice in the RG equations requires 
a new set of variables. Instead of the mass f, we use 
/o = 1/(1 -|- r) which corresponds to the on-shell prop- 
agator in the static limit, C(0), and asymptotically be- 
comes the inverse mass in the insulating phases where 
f(£) — >■ oo. In addition, we work with Ti — I^^i [i — 1, 2) 
for the temporal gradient terms and H — IqIi, G — I^g 
for the interaction and disorder vertices, respectively. In 
the following, we identify the MI phase by the conditions 
that f{l) oo but Rr{i) 0. In the BG, Rrii) ^ 0. 
In the large mass regime, G corresponds precisely to the 
relative variance Rr of the mass distribution. 



V. RESULTS 
A. Numerical Integration 

We start with a numerical integration of the full set of 
RG equations fTU]) in terms of the variables introduced 
above to see if we can indeed identify the MI/BG transi- 
tion. Since 71 7^ away from the tips of the Mott lobes, 
we use a dynamical exponent z = 2. Note that for small 
f close to the MI-SF transition in the clean system, the 
spatial gradient term 72; = 1 serves as a cut-off in C(lj) 
and the corresponding frequency integrals. 

Indeed we find a regime where the system flows to- 
wards a stable MI fixed point Pmi with Iq = and G = 0. 
Therefore, in the MI phase the mean f of the mass dis- 
tribution diverges fast enough that the relative variance 
Rr of the distribution renormalizes to zero. As expected, 
ri_2 and H renormalize to zero in the MI. 

Using the analytic expressions for the bare coefficients 
as a function of the microscopic parameters, we obtain 
the MI phases in the presence of disorder. Fig. 1(a) 
shows the resulting m — 1 MI lobe in = 3 as a function 
of x = ii/U and y = 2Dt/U for p = 0.1 and 6 = 0.2. 
Sufficiently close to the MI state, G diverges at a certain 
scale £* whereas Io{(*) is practically zero. This clearly 
identifies the adjacent phase as an insulating BG! Note 
that the mass distribution spreads faster than it shifts 
to infinity. Hence it is the tail of the distribution which 
destabilizes the MI, indicative of rare-region physics. 

In the following, we determine the phase boundary 
Ucix) by bisection and follow the divergence of the cor- 
relation length ^ = as y . Note that the corre- 
lation length defined here is determined by the scale £* 
where i?,- diverges and self-averaging breaks down. The 
correlation length ^ therefore corresponds to the distance 
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FIG. 1: (a) Phase boundary between the MI (m = 1) and 
the BG obtained from numerical integration of the RG equa- 
tions. The disorder corresponds to site energies increased or 
decreased by 5 = l\/U = 0.2 with probability p = 0.1. The 
solid line shows the unphysical mean-field phase boundary 
(f = 0) between the MI and SF. (b) Extraction of the corre- 
lation length exponent for different values x of the chemical 
potential as indicated in (a). 



of rare regions in the BG which cause the breakdown of 
self-averaging (see Sec. |lTl). Sufficiently close to the MI, 
we expect a power-law divergence, f ^ (y — yc)~'^ , defin- 
ing the correlation-length exponent v. To extract z^, in 
Fig. 1(b) we plot 9(ln^)/9(ln(y — j/c)) as a function of 
ln(j/ — yc) for different values of x. For y ^ y'^ , the func- 
tions indeed converge to a constant value corresponding 
to = 1/3. 

Interestingly, on approaching commensurate boson fill- 
ings, the length scale upon which universal behavior is 
observed increases and possibly even diverges. Whereas 
an increase might be simply a quantitative effect originat- 
ing from the suppression of the bare disorder vertex by 
commensuration, a divergence would indicate a different 
universality class at the tips of the MI lobes. 



B. Universality of the MI/BG Transition at 
Incommensurate Filling 

To better understand the nature of the transitions we 
proceed with an analytical investigation of the incom- 
mensurate case {z = 2). Since for D > 2 the interaction 
vertex h is irrelevant at the MI-SF transition of the clean 
system, in the vicinity of the fixed point Pmi-SF {G = 0, 
/q = 1) we can simplify the RG equations by setting 
h — 0. The resulting RG equations for Iq and G, 



dio 
d£ 
dG 

d7 



(G-2)/o + 2/2, 
{4Io-D)G + 6G^, 



(11a) 
(lib) 



are decoupled from Fi ^. Our analysis of these equations 
is strictly in the perturbative regime where G{£) < 1. 
The divergence of Iq which obtains for G{£) > 2 is spu- 
rious and outside the range of validity of the perturba- 
tive RG equations. In addition to Pmi-sf and the stable 
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FIG. 2: RG flow for (a) D = 3 and (b) D = 5 for incommen- 
surate boson fillings as a function of lo and G corresponding 
(in the limit of large mass f) to the inverse mean and the 
relative variance of the mass distribution, respectively. 



MI fixed point Pmi(^o = 0,6* = 0), the RG equations 
([TT1) exhibit an additional fixed point Pmi-bg(0, U/6), 
which will turn out to be critical fixed point determining 
the universality of the MI-BG transition at incommensu- 
rate fillings. From linearization around Pmi-bg we find 
dH/d£ — (-D/3 — 2)H indicating that the interaction ver- 
tex remains irrelevant. 

For G = we obtain a trivial fiow equation for Iq 
capturing the MI-SF of the clean system. For Iq < 1, the 
fiow is towards Pmi, while for /q > 1 (f < 0), /o(^) oo 
signaling the instability toward the formation of a SF. 

For D < 4, Pmi-sf is unstable for small disorder, 
G > 0. However, instead of a runaway fiow, we find 
a separatrix Gs(/o) which connects Pmi-sf and Pmi-bg 
as shown in the case D = 3 in Fig. 2(a). Lineariz- 
ing around the fixed-points, we find that asymptotically 
Gs(/o) ~ {D — 2)(1 — Iq) in the vicinity of Pmi-sf and 
Gs(/o) ~ D/Q ~ 4/o/(12 - D) close to Pmi-bg, respec- 
tively. 

For initial values /o(0) < 1 and G(0) < Gs(/o(0)), 
the system fiows towards the MI fixed point, whereas 
for G(0) > Gs(/o(0)), G(£) will eventually diverge. For 
initial values very close to the MI phase, the trajec- 
tory (/o(£),G(£)) will track the separatrix and conse- 
quently, sufficiently close to the transition, the divergence 
of G{1) is controlled by the fixed point Pmi-bg- Lineariz- 
ing around this point, we find G{t} = (AG)e^^ with 
AG — Go — Gs the infinitesimal distance to the sep- 
aratrix. From the condition G{£*) ~ 1 we obtain the 
correlation length f (AG)^^/-^ , implying a correlation 
length exponent oi v = l/D, in perfect agreement with 
the value extracted from numerical integration of the full 
set of RG equations for D = 3 (see Fig. 1(b)) and a 
previous real-space RG and scaling analysis. 

Since G 7^ at the fixed point and in the BG, the fiuc- 
tuations are no longer— governed by the central-limit 
theorem. This implies that self-averaging breaks down. 
Consequently, the bound 1/ >2/D [Ref. [2J| no longer ap- 
plies. That the BG is mediated by rare localized regions 
is the efficient cause of this effect. In quantum systems. 



such rare events play a more pronounced^S. role than in 
classical systems because of the perfect correlation of the 
disorder along the imaginary time direction. 

It is also the divergence of the relative variance 
that is responsible for a finite compressibility n = 
(/3C/)~^9^1nZ/9f^ in the BG. Using the asymptotic form 
of iScff in which only the r and g terms are retained, we 
obtain to leading order k w l/f^(l-|-4g/f^) which clearly 
vanishes in the MI. However, in the BG sufficiently close 
to the transition, k = k{1*) ^ ^0-4 ^ (^y _ with 
7 = A/D - 1. 

For completeness, we briefly discuss the case D > 4 in 
which a small amount of disorder is irrelevant at the MI- 
SF transition. The inversion of the RG flow along the G- 
direction is a consequence of the presence of an additional 
unstable flxed point Q{1 — e/8, e/4) in Z? = 4 -I- e located 
on the separatrix (see Fig. 2(b)). Therefore, in > 4 
the disorder has to exceed a critical value in order to 
induce a BG phase in contrast to recent claims 

C. Commensurate Filling 

We now focus on commensurate boson fillings to un- 
derstand the nature of the transition at the tips of the 
MI lobes where 71 = and z = 1. The smaller dynami- 
cal exponent makes the interaction vertex H — I^h more 
relevant and in fact marginal at the MI-SF transition of 
the clean system in D = 3. However, for D = 4 — e, 
H remains irrelevant at Pmi-sf suggesting that the crit- 
ical behavior is controlled by the same separatrix as for 
incommensurate fillings. Interestingly, this is not the 
case. A numerical integration of the RG equations shows 
that in the disordered phase close to the MI, H{1) di- 
verges simultaneously with G(Z). This behavior is eas- 
ily understood upon linearizing around Pmi-bg yielding 
dH/dl — {D/3 — 1)H which demonstrates that H be- 
comes relevant for Z) > 3. Therefore, for any finite value 
iZ(0) > 0, the separatrix obtained by projection into the 
/q-G plane does not terminate in Pmi-bg but diverges as 
Gs ^ Ig^^^'^^"^'. Consequently, the MI-BG transition 
at commensurate fillings is not controlled by Pmi-bg but 
by a different fixed point which is not accessible in the 
present calculation. However, from the divergence along 
the seperatrix, Gs ^ e"' we estimate v — Q/D which is 
significantly larger than the incommensurate value and 
satisfies the bound v > 2/ Di^ 

VI. FINAL REMARKS 

To summarize, we have analyzed the instabilities of the 
MI state of the weakly disordered BH model within an 
RG analysis of the disorder-averaged replica theory valid 
in the strong-coupling regime. Our analysis shows that 
the MI always becomes unstable towards a BG regard- 
less of the boson filling, in agreement with recent work— 
which, based on a mathematical theorem, excludes the 
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possibility of a direct MI/SF transition in the presence 
of disorder. 

While in both insulating phases the mass coefficient 
diverges on large scales, rendering SF correlations short 
ranged, the relative variance of the induced random mass 
distribution renormalizes to zero in the MI and diverges 
in the BG. This allowed us to demonstrate that the 
relative variance serves as the order parameter for the 
MI/BG transition. We further extracted the diverg- 
ing correlation length of the MI/BG transition from the 
scale where the relative variance diverges under the RG. 
Since the finite relative variance of an extensive quan- 
tity signals a breakdown of the central-limit theorem)^ 
the correlation length corresponds to the separation be- 
tween rare regions in the BG, causing a breakdown of 
self-averaging in the system. We point out that the cor- 
relation length defined here has nothing to do with the SF 
correlation length which remains finite over the MI/BG 
transition. 

The reformulation of the RG equations in terms of new 
variables, in particular the relative variance, allows us to 
identify fixed points and to subsequently analyze the uni- 
versality of the MI/BG transition. That new variables, 
not apparent in the UV-complete theory, are necessary to 
access the fixed point is typical of strong-coupling prob- 
lems. At incommensurate fillings, we find a separatrix 
which connects the MI/SF fixed point in the clean sys- 
tem with the critical MI/BG fixed point in the weakly 
disordered system. The correlation-length and compress- 
ibility exponents characterizing the MI/BG transition are 
determined as = l/D and 7 = 4/D — 1, respectively. In 
principle, the power-law divergence of the inverse com- 
pressibility could be verified either experimentally or nu- 
merically provided that the critical region can be ac- 



cessed sufficiently. The violation of the strict lower bound 
V > 2/D known as the quantum Harris criteriumSi has 
been argued some time ago^^^ to be a key signature 
for the breakdown of self averaging. Interestingly, our 
correlation- length exponent ly — l/D is identical to the 
exponent estimated by a scaling analysis in the atomic 
limit for a box distribution of site energies.-'^ 

Our analysis shows that the MI/BG transition at com- 
mensurate fillings is controlled by a different fixed point 
which is not accessible within the present RG approach. 
The two distinct transitions found here are, in retrospect, 
not out of the ordinary, since in the clean system, the 
MI/SF transitions are in the universality class of the 
{D + l)-dimensional XY modeP^^^* and of mean-field 
type- for commensurate and incommensurate fillings, re- 
spectively. This dichotomy persists in the disordered sys- 
tem as well. 

What this analysis indicates is that within a replica 
theory, instabilities driven by local rare regions, indica- 
tive of Griffiths singularities, are completely accessible 
if one focuses on the scaling dependence of the relative 
variance of the mass distribution. Since the relative vari- 
ance of any extensive quantity acquires a non-zero value- 
once self-averaging breaks down, the method presented 
here can in principle be applied to any transition driven 
by local rare-region physics. 
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